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Optical control of atomic interactions in a quantum gas is a long-sought goal of cold atom research. 

Previous experiments have been hindered by short lifetimes and parasitic deformation of the trap 
potential. Here, we develop and implement a generic scheme for optical control of Feshbach resonance 
in quantum gases, which yields long condensate lifetimes sufficient to study equilibrium and non¬ 
equilibrium physics with negligible parasitic dipole force. We show that fast and local control of 
interactions leads to intriguing quantum dynamics in new regimes, highlighted by the formation of 
van der Waals molecules and partial collapse of a Bose condensate. 


Spatiotemporal control of interactions would bring a 
plethora of new quantum-mechanical phenomena into the 
realm of ultracold atom research. Temporal modulation 
of interactions is theoretically proposed as a route for cre¬ 
ating anyonic statistics in optical lattices mm as well as 
new types of quantum liquids mm and excitations mm- 
Spatial modulation would grant access to unusual soli- 
ton behavior HE], controlled interfaces between quan¬ 
tum phases m, stable nonlinear Bloch oscillations El, 
and even the dynamics of acoustic black holes [12]. The 
conventional technique for controlling interactions in cold 
atoms, magnetic Feshbach resonance [13 HD, is typically 
insufficient for these applications because the magnetic 
coils are generally too large for very fast or local modu¬ 
lation. 

A promising alternative is optical control of Feshbach 
resonance (OFR). With laser beams, high spatial resolu¬ 
tion and high speed control of interactions can be realized 
by optical modulators. Efforts toward achieving OFR in 
quantum gases have made significant progress [T5H26] but 
encountered two major obstacles. First, in previous ex¬ 
periments OFR has limited the quantum gas lifetime to 
the millisecond timescale [20l |26] due to optical excita¬ 
tion to molecular states. Short lifetimes forbid studies of 
quantum gases in equilibrium or after typical dynamical 
timescales. Second, the change of interaction strength 
from OFR is often accompanied by an optical potential. 
This potential can result in a parasitic dipole force which 
dominates the dynamics when the interactions are spa¬ 
tially modulated [23] . 

In this report we propose and implement a novel 
scheme for optically controlling interactions while main¬ 
taining long quantum gas lifetime and zero parasitic 
dipole force. With a far detuned laser, a change of 
the scattering length a, which determines the interaction 
strength, by 180 Bohr radii (ao) is only coupled with a 
slow radiative loss of (0.63 s)“^. This loss rate is two to 
four orders of magnitude lower than previous reports with 
a similar change in scattering length and allows the BEC 
to remain in equilibrium. Eurthermore, the laser operates 
at a magic wavelength to eliminate the atomic dipole po¬ 
tential. We apply this OER scheme to test the response 
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FIG. 1. Illustration of optical control of Feshbach res¬ 
onances. A Feshbach resonance occurs when a laser (yellow) 
brings a molecular energy level (blue surface) close to the 
atomic scattering threshold (red surface). Here, the atom- 
molecule coupling makes atomic interactions more attractive 
at higher laser intensity. Operating at the magic wavelength, 
the beam does not shift the energy of single atoms (see text). 

of BECs to rapid oscillation of interactions down to the 
timescale of 10 ns, reaching beyond the van der Waals 
energy scale. By spatially modulating the interactions 
we observe intriguing dynamics, including the formation 
and collapse of solitons within a condensate. 

We optically control Eeshbach resonances by using a 
far-detuned laser to light shift molecular states near the 
atomic scattering threshold (EIG. [^. The large detun¬ 
ing from all atomic and molecular transitions offers low 
heating and loss rates for the quantum gas. Eor a laser 
with intensity /, the total light shifts of atoms (subscript 
a) and molecules (subscript m) are given by (see supple¬ 
ment) 

— (^a “1“ /^aMa)-^ 

where a is the scalar polarizability and the vector polar¬ 
izability Pii depends on the magnetic moment ja. Since 
our target molecular states are very weakly bound, they 
have similar polarizability to free atoms: a^n ~ 2(aa and 
Pm ^ /3a = /3 (supplement). Assuming the molecular 
and atomic magnetic moments differ fim 7^ 2/ia, the vec- 
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FIG. 2. Stable optical control of scattering length at a magic wavelength Am* a, Theoretical polarizability of 
Cs atoms in the absolute ground state for polarization (supplement). The red ic marks the magic wavelength where 
polarizability is zero, b, Measured polarizability (green o). A linear fit yields Am = 869.73(2) nm (red -k). c, Number of 
condensed atoms remaining over time with (red o) and without (blue □) exposure to OFR at a magnetic held of 48.19 G. We 
ht the decay dynamics (red and blue curves) and hnd that OFR exposure adds a one-body loss process with a time constant of 
0.63(2) s (supplement) at intensity I = 225 W/cm^. d, Scattering length a determined from the free expansion of BEGs with 
(red o) and without (blue □) exposure to the OFR laser. When the scattering length becomes negative (a < 0) the condensate 
collapses (x). The red and blue curves derive from a single ht to all a > 0 data using Eq. which yields A = 157(3) mG, 
Bo = 47.766(4) G, and j3I = —38(1) mG. e, In situ images of BEGs at 47.97 G (dashed line in d) after ramping on the OER 
intensity over 200 ms. Each image is the average of 10 trials. All error bars show standard error. 


tor light shift can bring the molecular states closer to the 
scattering state, inducing a resonant atom-molecule cou¬ 
pling. Moreover, we choose a magic wavelength Am to 
eliminate the dipole force on the atoms = 0) [27] . 
such that only the molecular shift 

« ^(/i^ - 2/ia)/ (1) 

remains (FIG. [^. Under these conditions the laser can 
change the scattering length without creating parasitic 
dipole forces. This scheme can be implemented in atomic 
species with a magnetic Feshbach resonance and a magic 
wavelength far detuned from electronic transitions (sup¬ 
plement). 

We implement this OFR scheme in cesium BEGs pre¬ 
pared with a small positive scattering length near the 
Feshbach resonance at 47.8 G (see supplement for exper¬ 
iment details). Theoretically, one possible magic wave¬ 
length is 869.7 nm for a polarized laser, which is 
far detuned from all electronic transitions but main¬ 
tains a large vector polarizability (FIG. If). We con¬ 
firm the magic wavelength experimentally by measuring 
the dipole force of the OFR laser on the atoms (sup¬ 
plement). From the extracted polarizability, we deter¬ 
mine the magic wavelength to be Am = 869.73(2) nm 
(FIG. [^). At the intensity / = 225 W/cm^ used for 
most of this work, we estimate that the residual dipole 


potential /cp x 1 nK is negligible compared to our typical 
chemical potential of /cp x 10 nK, where /cp is the Boltz¬ 
mann constant. With uniform exposure to this intensity, 
the loss induced by OFR is well explained by a one-body 
time constant of 0.63(2) s (FIG.[^). 

To precisely determine the change of scattering length, 
we allow the BEG to freely expand with and without 
exposure to the OFR laser and measure the size of the 
gas after expansion [28l[29] (supplement). FIG.[^ shows 
the shift of scattering length induced by the laser near 
the 47.8 G Feshbach resonance. We fit the scattering 
lengths with a theoretical model na 


(2(/) — (2bg 



A 

B[I) - 5o 


( 2 ) 


where B{I) = + (31 is the effective magnetic field 

including the OFR contribution (Eq. Q and the ex¬ 
ternal field Bex, CL]jg ^ 950 ao is the background scat¬ 
tering length at this Feshbach resonance [30], A is the 
width of the resonance, and Bq is the resonance position. 
The fit yields pi = —38(1) mG, sufficient to decrease the 
scattering length from 180 ao to zero. 

The long BEG lifetime allows us to corroborate the 
change of scattering length based on in situ measure¬ 
ments of the density profile. We slowly ramp on the OFR 
beam to four different final intensities over 200 ms and 
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FIG. 3. Interaction modulation spectroscopy. The BEG at 47.976 G is exposed for a time t to the OFR beam, 
which is intensity modulated at frequency ujI'I'k. The number of condensed atoms remaining after exposure, normalized to 
the off-resonant number, is measured under three conditions: (o) t — 100 ms with no optical lattice; (•) t = 20 ms with a 
one dimensional optical lattice of depth h x 9.3 kHz; (inset □) t = 500 ms with no optical lattice. Resonances are observed 
at 0.89(1), 8.18(2), 133(7), 197(1), and 11,649(2) kHz, determined from hts (solid curve) to each resonance (Gaussian for 
133 kHz, Lorentzian for others). The illustrations indicate the nature of each resonance. The van der Waals energy scale is 
E^vdw — h x2.7 MHz for Gs molecules. 


measure the resulting column density profiles (FIG.[^). 
Higher OFR intensities shrink the BEG and increase its 
density, consistent with weakening the repulsive interac¬ 
tions. A mean-field analysis yields scattering lengths in 
excellent agreement with the free expansion measurement 
(supplement). 

The stability of this scheme enables us to explore tem¬ 
poral and spatial control of interactions in a quantum gas. 
We first perform interaction modulation spectroscopy by 
recording the response of the BEG to an OFR beam with 
oscillating intensity, see FIG. We observe a variety of 
resonance features over a wide range of timescales, high¬ 
lighting the versatility of this technique. 

At lower frequencies we observe excitations in the trap 
and in the optical lattice which directly result from the 
oscillating scattering length. This demonstrates that 
plentiful theoretical proposals which require rapid 
oscillation of scattering length in the bulk or in the lattice 
are well within reach of our scheme. 

At higher frequencies the oscillating OFR intensity 
induces formation of molecules. We identify reso¬ 
nances corresponding to a virtual state at 133(7) kHz, 
a weakly-bound Feshbach molecular state at 197(1) kHz, 
and a deeply-bound van der Waals molecular state at 
11.649(2) MHz. All observed resonance positions are in 
excellent agreement with theoretical calculations (sup¬ 
plement). These resonances provide direct evidence that 
OFR can access interaction physics on timescales as short 
as 10 ns. 

Next, we demonstrate spatial modulation of the inter¬ 
action strength within a quantum gas. For this experi¬ 


ment we employ an OFR beam which is small compared 
to the size of the BEG, leading to a reduced scattering 
length only in the center of the gas. After preparing the 
BEG we quickly turn on the OFR beam and simultane¬ 
ously change the magnetic field. We study the subse¬ 
quent dynamics of the sample by measuring its in situ 
density profile over time (FIG. |^. For example, when 
the interactions remain repulsive throughout the conden¬ 
sate (FIG. 1^), we observe collective excitations for the 
duration of the experiment. For all images the small dis¬ 
tortion at the center of the gas results from the dipole 
potential due to slightly non-uniform laser polarization 
(supplement). 

Intriguing quantum dynamics occur when the inter¬ 
actions become locally attractive. When the scattering 
length is negative in a small region near the center of the 
trap (FIG. [^), we observe a brief period of transverse 
compression followed by a rapid drop in central density 
between 20 and 32 ms after the quench, signalling lo¬ 
cal collapse of the condensate (illustrated in FIG.[^). A 
large fraction of the sample survives at the edges for more 
than 100 ms. With even stronger attractive interactions 
(fig.| 5), faster central collapse occurs after 8 ms. Sub¬ 
sequently, the dense remnants at the edge of the sample 
move toward the center of the trap (see white dashed 
lines in FIG. [^). Beyond 32 ms only thermal gas sur¬ 
vives, indicating that the remnants have undergone fur¬ 
ther collapse. Based on this behavior, we identify these 
remnants as bright matter wave solitons [31] , which form 
at small negative scattering length but become unstable 
when the scattering length drops below a critical value. 
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FIG. 4. Condensate dynamics with spatially mod¬ 
ulated interactions. Time series of in situ images of 
BECs with N = 12000 atoms after a quench from uniform 
a = 200 ao at 47.965 G to the spatially modulated a shown 
in the top panels. The OFR beam has peak intensity of 
115 W/cm^ and a waist of 14 /xm, while the final magnetic 
fields are a, 47.949 G, b, 47.935 G, and c, 47.925 G. Each im¬ 
age is the average of 6 or 7 trials. The red dashed lines show 
where a equals zero. The white dashed lines in c guide the eye 
to the motion of the solitonic wave towards the trap center, 
d, Illustration of the local collapse dynamics, in which the 
initial BEG (blue) undergoes transverse compression followed 
by localized central collapse (yellow). 


The variety of behaviors observed in this experiment es¬ 
tablishes the richness of the quantum dynamics accessible 
with space dependent interactions. 

In conclusion, we implement a generic scheme for op¬ 
tically controlling interactions in quantum gases. This 
scheme overcomes the key challenges encountered by past 
approaches to OFR. Fast and local control of interac¬ 
tions in a quantum gas enables studies of novel quantum 
dynamics and has great potential in the fields of quan¬ 
tum simulation and engineering, opening a new frontier 
of interaction-driven quantum phenomena. 
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SUPPLEMENTARY INFORMATION 
Experiment Setup 

Our experiments test optical control of Feshbach res¬ 
onances in almost pure BECs of cesium atoms formed 
in a crossed beam optical dipole trap (wavelength A = 
1064 nm). This harmonic trap has typical horizontal 
frequencies (uJx^oOy) = 27r x (12,30) Hz and vertical fre¬ 
quency ujz controlled by the beam intensities, which we 
vary from 27r x 70 Hz to 27r x 470 Hz to suit the needs 
of individual experiments. Our samples typically con¬ 
tain from 3 r\j 10 xlO^ atoms with peak densities from 
1^4x10^^ cm“^. An objective lens (NA=0.5) oriented 
vertically collects our in situ absorption images with 


^1 fim resolution. When we test the response of con¬ 
densates in optical lattices to oscillating OFR intensity, 
one of the horizontal trapping beams is retro-reflected to 
create a ID optical lattice with spacing d = 532 nm and 
depth h X 9.28 kHz. For most of the experiments de¬ 
scribed in this paper we form our BECs near the d-wave 
magnetic Feshbach resonance at 47.8 G [32] with a small 
and positive scattering length of a = 200^300 uq. The 
magnetic field in our system is stable to within 1 mG 
with systematic calibration error of 5 mG, sufficient to 
form stable BECs. 

We optically control interactions with the intensity of 
the OFR laser, which propagates along the magnetic field 
direction. The change of scattering length due to OFR 
is insensitive to small changes in the wavelength and po¬ 
larization of the OFR laser, however both parameters 
must be carefully controlled to eliminate the residual 
light shift. The source of the laser is a free-running single¬ 
mode diode with 100 mW output power, which can be 
tuned via its temperature. We determine the laser fre¬ 
quency using a wavemeter with an accuracy of 100 MHz, 
sufficient to make the residual light shift negligible. 

The temporal modulation experiment (see Fig. 3 in 
main text) requires rapid control of the beam intensity. 
We send the OFR beam through an acousto-optic mod¬ 
ulator (Isomet 1205C-1) and use a fast RE switch (Mini- 
Circuits ZFSWA-2-46) to control the acoustic wave which 
drives the modulator. We measure the 3 dB bandwidth 
of intensity modulation to be 10 MHz. 


Theoretical calculation of light shift, scattering rate 
and performance of OFR 

To predict the magic wavelength Am and the ef¬ 
fective field shift pi from OFR (main text Eq. 1) 
we calculate the scalar and vector polarizabilities. 
Our Cs atoms are prepared in the absolute hyper- 
fine ground state |F = 3, Mp = 3) where the total an¬ 
gular momentum F = J + / is the sum of elec¬ 
tron angular momentum J = 1/2 and nuclear spin 
/ = 7/2, and Mk {K = F, J) is the projection 
onto the quantization axis. At low field we have 
\F = 3, Mf= 3) = -^/7j8\Mj = -1/2, Mi = 7/2) + 
v/T78|Mj = 1/2, Mi = 5/2). 

The scalar (a) and vector {/dja) AC polarizabilities of 
an atom for detuning which is large compared to the 
hyperfine splitting are [33] 
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Wavelength (nm) 


its expectation value of < Mj >= —3/8 for atoms in the 
|F = 3, Mf = 3) state. The calculated polarizabilities 
are shown in FIG. [STl 

The figure of merit M for choosing the wavelength and 
polarization of the laser is the ratio of the effective field 
shift pi from OFR to the photon scattering rate s{I)^ 



8{I) = 


I sr 

2ffieoc ^ {uifi - ujY ’ 


(S3) 


FIG. SI. Theoretical polarizability. Scalar (black), (t"*" 
vector (red dashed), and cf~ vector (blue dot-dashed) polar¬ 
izabilities in the absolute ground state of Cs, see equations 


are included in the calculation. Two possible magic wave¬ 
lengths exist with circularly polarized light, 869.7 nm for cr^ 
and 891.0 nm for cf~ polarization. We employ polarization 
for this work. 


SI and |S2[ Only the D1 (894 nm) and D2 (852 nm) lines 
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where |/) represents a relevant excited state, eo is the vac- 

I 71 72 7s 

uum permittivity, c is the speed of light, < . . . 

t 3b J6 

is the Wigner 6-j symbol, |(/||d||i)| and ujfi are the re¬ 
duced dipole matrix element and resonance frequency for 
the transition from \i) to |/) (Table [^, uj is the laser fre¬ 
quency, Mj is the projection of Ji onto the laser prop¬ 
agation direction, and A = (I^+ — I(j-)jl accounts for 
the beam polarization. We include in the calculation 
only the D1 and D2 transitions, which dominate in the 
wavelength range that we consider. The next most sig¬ 
nificant excited state provides a correction at the 10“^ 
level, which is negligible for our purposes. We set Mj to 


where = (7^/(/||d||i)p is the squared dipole matrix 
element, (7^^ is a numerical factor primarily accounting 
for the Clebsch-Gordan coefficients with beam polariza¬ 
tion , and Fj is the spontaneous emission rate of the 
excited state |/) (Table 0. A larger absolute value of M 
indicates a greater shift in molecular states for a fixed 
quantum gas lifetime. 

M is optimized by using pure circular polarization. 
The magic wavelengths for and a~ polarization are 
both suitable choices with almost identical M to within 
5%. They remain differentiated by the direction of the 
change in scattering length, dajdl < 0 for cr+ and 
da/dl > 0 for a~. We employ cr+ polarization with M = 
— 130 mG-s at the magic wavelength Am = 869.7 nm. 

Optimal performance of our OFR scheme is obtained 
when a laser with maximum M is employed near a proper 
magnetic Feshbach resonance. The sensitivity of scatter¬ 
ing length to the OFR shift (based on Eq. 2 of the main 
text) is 


da 


^bg 


I Q'/O'bg I '^1 


A 


(S4) 


near the zero crossing B{I) = Bq-\- A where the conden¬ 
sate can be stable. This result suggests that a narrow 
resonance with a large background scattering length will 
offer wide tunability for a given laser intensity. However, 
one must be cautious to avoid extremely narrow reso¬ 
nances which can greatly enhance the three-body recom¬ 
bination rate. In this work we choose the Feshbach reso¬ 
nance at 47.8 G, for which we have Ubg/A = 6.6 ao/mG 
from coupled channel calculations m- 

Overall, optimal performance is obtained via the inde¬ 
pendent optimizations of M, which is determined by the 
laser, and dajdi^pip which is determined by the choice of 
Feshbach resonance. In our system the product of these 
two factors yields M(abg/A) =—860 ao-s. This value 
indicates that a decrease in scattering length by 860 ao 
should be possible with a scattering rate of 1/s. 
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ujfil27T (THz) 

|(/||d||*)| (10-2® C-m) 

Tf /2Tr (MHz) 

F 


D1 

335.116 

3.81 

4.557 

7 

24 

1 

24 

D2 

351.726 

5.36 

5.219 

5 

48 

11 

48 


TABLE L 


The parameters used to calculate the polarizability and scattering rate for the D1 and D2 transitions of Cs. 


Implementing OFR with other atomic species 

Our scheme for OFR is quite general and should per¬ 
form well with a variety of atomic species. For example, 
calculations based on ^^Rb in the absolute ground state 
\F = 1, Mp = 1) yield the figure of merit M = —45 mG-s. 
Moreover, past experiments with ^^Rb suggest two po¬ 
tential candidates for the Feshbach resonance. For ^^Rb 
in \F = l^Mp = 1), the 1007.4 G resonance [19] has a 
sensitivity of dajdifiT) = 0.5 ao/mG [T4|, which would 
allow a change of scattering length by 100 uq with a 
5/s one-body scattering rate. If a magic wavelength is 
not required, then the interstate Feshbach resonance at 
9.13 G |34| offers dajd^^F) — 6.7 ao/mG, where a is 
the scattering length between the |F = 1, Mp = 1) and 
|F = 2, Mp = — 1). This resonance allows a change of 
300 ao with a 1 s“^ scattering rate, enabling many inter¬ 
esting studies with optical control of interstate scattering. 

Beyond alkali atoms, quantum gases of heavy, mag¬ 
netic atomic species like erbium and dysprosium [35j [36| 
are promising candidates for OFR. Their rich sets of op¬ 
tical transitions provide many options for obtaining a 
favorable ratio of the vector light shift to the scatter¬ 
ing rate. Moreover, such highly magnetic elements have 
an abundance of Feshbach resonances naiiH], many of 
which could provide favorable properties for implement¬ 
ing OFR. 


Measurement of polarizability 


We find the magic wavelength to be Am = 
869.73(2) nm for polarization, in fair agreement with 
the prediction of 869.66 nm. We attribute the difference 
to the electronic transitions not included in the calcula¬ 
tion and the imperfect beam polarization; we estimate 
that Icr- « 0.005/. 

In addition, at the magic wavelength we test the as¬ 
sumption that 13m ^ Pa = P for weakly-bound molecular 
states. We measure the effective field shift pi induced by 
the laser using microwave spectroscopy. We compare the 
microwave result to the shifts in energy of g-wme and 
d-wave molecular states. Within 5% all methods yield 
values which are consistent with the theoretical predic¬ 
tion (Eq. [S^. 

Extracting scattering length from BECs after free 
expansion 

A precise way to measure the scattering length is to 
release the BEG from the harmonic trap, simultaneously 
switch the magnetic field and OFR beam intensity, and 
measure the BEG radius after a period of free expan¬ 
sion. This method yields measured radii R which scale as 
a^/^, making it more sensitive than in situ measurements 
or free expansion measurements with constant magnetic 
field which yield R oc a^/^ [29]. Ref. [28] shows that, dur¬ 
ing expansion, the BEG density profile is parabolic with 
time-dependent Thomas-Fermi radii Rk {k = x^y^z). 
The time evolution of the ratios \k{t) = Rk{t)/Rk{0) 
is described by the three coupled differential equations. 


We determine the polarizability near the magic wave¬ 
length by measuring the displacement of the BEG caused 
by the OFR laser. For this experiment we operate at 
B = 22 G where the dependence of interactions on the 
OFR shift is negligible. We prepare the OFR beam 
with a waist of 70 /im, which is large compared to 
the BEG, and displace the beam to achieve a slope of 
dl/dx = b (W/cm^) //im across the sample. A finite total 
polarizability j = aF Py leads to a force of F = jdljdx 
on the atoms, which offsets the center of the harmonic 
trap by Ax = F/muj‘^. By measuring the shift in the 
center of the BEG, we extract the polarizability 


mu)‘^Ax 
^ dl/dx ^ 


(S5) 


yielding the data shown in FIG. 2b of the main text. 


A/c 




(S6) 


where ujk is the trap frequency before release, a^ = 250ao 
is the scattering length before release, and a/ is the scat¬ 
tering length during expansion which is determined by 
the magnetic field and OFR intensity. Note that there 
is one equation for each A/^ and the trap frequencies are 
independently calibrated, such that the only unknown 
variable is the scattering length. 

To determine a/ as a function of magnetic field and 
OFR intensity, we execute the experimental procedure 
over a range of magnetic fields with and without ex¬ 
posure to the OFR beam at / = 225 W/cm^. For 
this experiment we use a vertical trapping frequency of 
ujz = 27 r X 105 Hz much greater than the horizontal fre¬ 
quencies of (ujx^ooy) = 27 r X (14, 31) Hz, which causes the 













Magnetic Field - 48G (mG) 

FIG. S2. Comparison of scattering length from free 
expansion to coupled channel model. The scattering 
length extracted from free expansion with no OFR laser (blue 
o) compared to the scattering length based on the coupled 
channel calculations in Ref. [30] (solid curve). The dashed 
curve shows the calculated scattering length shifted by - 
16 mG, which is within the model uncertainty. 



Column Density (pm'2) Position (|jm) 


FIG. S3. Effect of OFR on in situ density profiles. 

a, The in situ image series from fig. 2e of the main text. 

b, Horizontal line cuts (•) obtained by averaging the middle 
five pixels from each image in a, similarly arranged from top 
(no OFR laser) to bottom (225 W/cm^). The solid lines are 
the corresponding cuts from a fit to each image with Eq. [S7| 
For the 225 W/cm^ case we first perform a Gaussian fit to the 
wing to account for thermal atoms before fitting the remaining 
BEG. The fitted scattering lengths, in the order of increasing 
beam intensity, are a = 230(33), 175(25), 120(16), 96(13), and 
78(16) ao in good agreement with the free expansion result. 


vertical expansion to dominate after release. Thus, after 
an expansion time of 16 ms we measure only the vertical 
radius of the gas. For each combination of magnetic field 
and OFR intensity, we determine the scattering length af 
by comparing the measured radius to the numerical so¬ 
lution of Eq. This procedure yields the results shown 
in FIG. 2d of the main text. In the absence of the OFR 
beam we compare our results to coupled channel calcu¬ 


lations [30|, see FIG. S 2 


Extraction of scattering length from in-situ BEG 
density profiles 

We fit the measured in-situ column density of BECs 
with a Thomas-Fermi density profile [39] 


/ 2 2 \ 

= , (S7) 

/ 2 \ 2/5 

where no = (^) 

column density, and Rk = ‘ 

The only free parameter is the scattering length a which 
appears in no and Rk^ while the atom number N = 3500 
and trap frequencies {ux^ojy^ujz) = 27r x (10,22,75) Hz 
are calibrated independently. We compare horizontal line 
cuts of the images to the corresponding cuts from the fit 
in FIG. ESI 


Condensate lifetime with OFR 

In the presence of the OFR laser, we expect the life¬ 
time to be limited by one-body off-resonant scattering 
of photons, since our scheme does not rely on proximity 
to resonant atomic or molecular transitions. In prac¬ 
tice, while we predict the one-body loss process to have 
extremely weak wavelength dependence, we find signif¬ 
icant wavelength dependence of the lifetime within the 
40 GHz tuning range of our laser. We attribute this ob¬ 
servation to the two-body loss caused by photoassocia¬ 
tion or molecular resonances which happen to be within 
the tested wavelength range. After finding a wavelength 
which minimizes this loss, we make minute adjustments 
to our laser polarization to shift the magic wavelength to 
the loss minimum. The required change in polarization 
is small, such that its effect on the effective field shift is 
negligible. 

Having minimized two-body loss, we model the lifetime 
by accounting for the one-body loss process induced by 
OFR and the pre-existing three-body loss process. With 
both processes the decay of density n is described by the 
differential equation 

=-Lin{r,t) - Lzn{r,tf, (S 8 ) 

where Li (I/ 3 ) is the one-(three-)body decay constant. 
We assume that the loss is slow such that the system 
remains in equilibrium. Integrating the Thomas-Fermi 
density profile [39] over space with fixed scattering length 
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FIG. S4. Energy levels for Cs 2 molecules. A schematic 
of the three molecular resonances observed in our experiment, 
based on the energy levels in ref. [32]. The red ball indicates 
the condensate, and arrows represent the coupling to three 
molecular states by interaction modulation spectroscopy. 


and absorbing trap-dependent factors into the modified 
decay constant L'^ we obtain, 


^ = -LiN - 
dt ^ 


(S9) 


which describes the decay of total atom number in the 
presence of one- and three-body loss. Note that the one- 
body lifetime is independent of trap geometry. 

To determine the lifetime experimentally, we compare 
the decay of condensate number with and without ex¬ 
posure to the OFR laser at = 48.19 G (main text 
FIG. 2 c). We chose this field to obtain realistic decay 
rates near the Feshbach resonance but be far enough 
away to avoid inducing condensate dynamics or signif¬ 
icantly altering the three-body loss with the OFR beam. 
We keep the trap depth constant throughout the decay, 
so that the BEC continues to evaporate and maintain 
an approximately constant temperature while exposed 
to the laser. With no OFR exposure we fit the decay to 


the numerical solution of Eq. S9 with negligible Li and 
find I /3 = 5.6 X 10 ^ s ^ in our trap geometry. We use 
this fixed value to fit the decay in the presence of OER 
to Eq. S9 and obtain Li = 1.58(4) s“^, corresponding 
to a time constant of r = l/Li = 0.63(2) s, with OER 
intensity 225 W/cm^. 

The fitted time constant represents two orders of mag¬ 
nitude improvement over existing OER schemes, but it 
is still shorter than the off-resonant scattering time con¬ 


stant of 3.2 s (Eq. S3). We expect that the lifetime is 
shorter than the scattering time constant because the re¬ 
coil heating from a single scattering event is x 200 nK, 
much larger than the critical temperature Tc ~ 30 nK 


of the condensate. Thus each scattering event can lead 
to the loss of multiple atoms from the condensate. There 
may also be a contribution from photoassociation reso¬ 
nances or molecular transitions which remain near the 
chosen wavelength. 


Theory comparisons for temporal modulation 
resonances 

We have determined the origin of each of the reso¬ 
nances excited by interaction modulation spectroscopy 
(main text EIG. 3). The first resonance at 890(10) Hz 
corresponds to excitation in the vertical harmonic trap. 
Since the BEC density is even across the trap center, the 
oscillating interaction strength provides an even pertur¬ 
bation which can only excite the gas to states with the 
same parity as the ground state; thus the first excited 
state is forbidden. The ratio of the resonance frequency 
to the vertical trapping frequency of 470 Hz is 1.9, sug¬ 
gesting that the samples are close to the quasi-2D regime 
(chemical potential fi 470 Hz) and OER is driving ex¬ 
citations to the 2 nd excited harmonic oscillator state. 

The second resonance at 8.18(2) kHz corresponds to 
excitations to the second excited band of the optical lat¬ 
tice. Eor comparison we perform band structure calcu¬ 
lations for our lattice depth of h x 9.28 kHz, yielding 
energies of x 6.58 kHz and h x 8.18 kHz for the zero 
quasimomentum states in the first and second excited 
bands, respectively. The observed resonance is consistent 
with transitions to the second excited band; no resonance 
appears at the first excited band energy because the par¬ 
ity of the ground state is again conserved. Transitions to 
higher parity-allowed bands are also not observed, which 
we attribute to the weaker coupling strength to those 
states. 

The highest three resonances correspond to coupling 
of the free atoms to Cs 2 molecular states. The process 
is analogous to the binding of dimers with an oscillat¬ 
ing magnetic field [40]. The locations of the observed 
features are consistent with known Cs 2 molecular states 
(EIG. [S4| ). We compare our data to coupled channel 
calculations m at the average effective magnetic field 
B = 47.957 G. The 133(7) kHz resonance corresponds to 
a molecular state above the atomic threshold, with an en¬ 
ergy of 95 kHz based on extrapolation of coupled channel 
calculations into the continuum. The feature is broad be¬ 
cause the molecular state is embedded in the continuum. 
The 197(1) kHz resonance corresponds to a primarily d- 
wave bound state calculated to be at 208 kHz and the 
11.649(2) MHz resonance corresponds to a ^-wave bound 
state calculated at 11.731 MHz [30]. These molecules are 
individually stable, but their collisions with each other 
and the remaining atoms lead to net loss of atoms from 
the trap. 












